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WITHOUT AVERAGING 
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Abstract. We prove the existence of a large class of one-parameter families of cosmological solutions 
to the Einstein-Euler equations that have a Newtonian limit. This class includes solutions that repre- 
sent a finite, but otherwise arbitrary, number of compact fluid bodies. These solutions provide exact 
cosmological models that admit Newtonian limits but, are not, either implicitly or explicitly, averaged. 



1. Introduction 
Gravitating perfect fluids are governed by the Einstein-Euler equations 

G ij = ^±T 3 and V;T« = 0. 
c 4 

where 

= {p + c- 2 p)v l v J +pg l \ 

with p the fluid density, p the fluid pressure, v l the fluid four-velocity normalized by v z Vi = — c 2 , c 
the speed of light, and G the Newtonian gravitational constant. By suitably rescaling (see [8]), these 
equations can be written as 

(1.1) G 13 = 2e 4 T i] and V,T 4J = 0, 
where 

T %3 1 = (p + e 2 p)v V ; + pg ij and v l Vi = —-^. 

In the article [10], we proved the existence of a large class of one-parameter families of solutions to this 
system defined for < e < eo that 

(i) exist on a common piece of spacetime of the form M = [0, T) x T 3 , 

(ii) converge as e \ to solutions of the cosmological Poisson-Euler equations of Newtonian gravity, 
and 

(iii) are differentiable in e to any prescribed order I 6 N. 

The properties (i)-(iii) guarantee that these one parameter families of solutions to the Einstein-Euler 
equations have valid Newtonian limits and admit post-Newtonian expansions to order £/2. However, 
in order to establish the existence of these solutions, we required that the two following conditions are 
satisfied by the initial fluid density p\t=o = po and fluid three- velocity u / |t=o = Wq'. 

(1.2) Po{x)>0 VieT 3 , 
and 

(1.3) / p w I d 3 x = Q. 
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Here, we are using x = (x 1 ) (I = 1,2,3) to denote the standard periodic coordinates with period 1 
on the torus T 3 = S* 1 x S 1 x S , and t = x° to denote an absolute Newtonian time coordinate on the 
interval [0,T). 

The main aim of this article is to remove the conditions (|1.2[) - (|1.3[) on the initial data. Although, we 
do not discuss the existence of post-Newtonian expansions in this article, it is not difficult to see that 
the results of this article can be combined with those of [10] to prove the existence of post-Newtonian 
expansions to arbitrary order that do not satisfy l|1.2|l - (|1.3|) . 

From our point of view, there are two important reasons for removing the conditions (|1.2[) - p.3p . This 
first reason is that (|1.3p is an exact averaging condition that cannot be expected to be exactly satisfied 
for real systems. The second is that by removing (|1.2p , we can choose an initial fluid density of the form 

JV 
u=0 

where N is any integer, and 

dist(supppo,i/,suppp ,p) > v 7^ M- 
This initial data represents a finite but otherwise arbitrary number of compact fluid bodies (i.e. stars). 
Thus, solutions generated by this initial data can be used to model an arbitrary collection of stars, and 
therefore provide an exact model for the universe. Consequently, we obtain one-parameter families of 
cosmological solutions that admit Newtonian limits, and are not, either implicitly or explicitly, averaged. 
To formulate the Newtonian limit, we require (see [10]) the following FLRW dust solution = 
C,p e =^Pe = 0) of (□}: 



(1.4) 
(1.5) 

(1.6) 



^ = 



a £ S(Sj8 



ij. 



where a e — a e (t) and [i e — /j, e (i) satisfies the differential equations 

/ 1 2 



(1.7) 

and 

(1.8) 

respectively. Here, we are using 



, 3 a' 
2 a f 



;Me = 0, 



<■»' - !«• 



The differential equations (|1.7[) - (|1 ,8f> can be integrated explicitly to give 
(1.9) 



a £ (t) = a e (0)(./^P < + i 



and 
(1.10) 
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where a e (0) and cr(e) = /i e (0) are arbitrary (positive) functions that are analytic in a neighborhood of 
e = 0. We fix the arbitrary length scale by setting a e (0) = 1, while a(e) will be determined by later 
considerations. Also, to reduce notation we will drop the e from the function a and \x except in situations 
where we want to emphasize the e-dependence. 

In this article, we use a slight variation of the approach used in [10] to analyze the limit e \ 0. The 
first step in the analysis is to replace the metric and the fluid velocity v l with variables that are 
compatible with the limit e \ 0. The new gravitational variable u lJ is defined by 

< U1 > s, ' = ^Jko> <w = 

where 

(1.12) q 11 = h ij + e 2 JlJ{u kl and j{ = eSffi + SjSf, 
while the new fluid four-velocity to' is defined by 

(1.13) w* = (1 + ew°)5i + Sl-w 1 . 

As in [8,9], we use Makino's technique [7,11] to generate perfect fluid solutions with compact support. 
This requires the use of an isentropic equation of state of the form 

(1.14) p = Kp {n+1)/n , 

where K £ M>o, n G N. This allows for the introduction of the density variable 

^' 15 ^ P ~ (4Kn(n+l)) nQ 

which is used to formulate the Euler equations as a symmetric hyperbolic system that is regular across the 
fluid- vacuum interface. In this way, it is possible to construct solutions to the Einstein-Euler equations 
that represent compact gravitating fluid bodies (i.e. stars) both in the Newtonian and relativistic 
setting [7,11]. Although, these solutions do not include static stars of finite radius [11], they are general 
enough to understand the mathematical issues involved in the Newtonian limit. 

The main result of this article is to show that solutions to the following equations, which we refer to 
as the cosmological Poisson-Euler-Makino equations, rigorously approximate fully relativistic solutions 
to the Einstein-Euler equations (jl.ip up to an error term of order e as measured in suitable Sobolev 
spaces: 

(1.16) dta — — w 1 dia — ^-diw 1 — -7— -a, 

2n Ana 

(1.17) d t w J = -^d^- 1 - -^r<9 J d - ^rw J + g J , 

2na a 

(1.18) A$ = 4505-/2) (p = (AKn(n+ 1))~"<S 2 "), 
where 

a = a , /2 = fi , 

(1.19) g J = — (--- f pw J d 3 x+-d J ^ 

a \ 2 a p, J T 3 4 

and (-|-)l 2 is t ne standard L 2 inner-product on T 3 , i.e. 

(^iIV^)^ = / ipi(x)tl) 2 (x)d s x. 

J[0,1] 3 

Here, and for the rest of the article, 

A = 5 IJ didj 
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will denote the flat Laplacian, and we set 



d J = S IJ dr. 



7- 

In terms of the fluid density p, equations (|1.16|) - l|1.18|) take the form 

3 a f 

(1-20) d t p = -w J dip- pdjw 1 - -—p, 

2 a 

(1.21) d t w J = -w z diw J - ^zd J p(p) - ^rw J + g J , 

ap a 

(1.22) A$ = 4a(/5-/2), 

which we refer to as the cosmological Poisson-Euler equations. We note that these equations agree with 
the Newton-Cartan field equations for a gravitating fluid formulated in adapted coordinates [5, 12]. 

An important point is that the constant er(0) = p(0) is not arbitrary but is determined by the initial 
density 

Pa = p\t=o 

according to 



(1.23) (To := A(0) = / Pod 3 x > 



<jo ■= £(o) = / 

Jt 3 



From equations (|1.8|) and (|1.20[) , it is not difficult to see that this averaging condition is preserved under 
evolution 



(1.24) p(t) = / p(t)d 6 x, 

which in turn implies (see (|1.6|) ) that 



(1.25) a(i)=exp / -/ p{s)d 6 x\ ds 



o V3 



T 3 



Also, setting 

(1.26) C J = / pw J d 3 x 

JT 3 

it follows directly from the cosmological Poisson-Euler equations (|1 .20[) - (| 1 . 22[) that ( J satisfies the equa- 
tion 



a 



which can be integrated to give 

(1-27) C'(t) = ^j f Po^d 3 x, 

a(t) J T 3 

where 

w 7 = w J \t=o 

is the initial fluid 3- velocity. Using Q1.27p . the acceleration due to gravity (|1.19[) can be reexpressed in 
the form 

a \a'a 6 J T 3 4 
For purposes of interpretation, it is often useful to introduce Galilei coordinates [5,12]. These co- 
ordinates are defined as follows: suppose {p(t, x), u/(i, x), <fr(t, x)} is a solution of the cosmological 
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Poisson-Euler equations fOOjl - fO^]) on M = [0, T) x T 3 . Then, letting M = [0,T) x R 3 denote the 
covering space, we define a diffeomorphism on M by 

■0 : M — >M : (t,ar) i — ► (i, x/s/d^j). 

Lifting the cosmological Poisson-Euler equations to M, and then pulling back by ip shows thalQ 

(1.28) p(t,x)=p(t,x/y/di(j$), 

(1.29) w J (t,x) = % /^w J (t,x/^/U£)) + \—^- xJ ^ 

(1.30) Ht,x) = §(t,x/^/5tfj), 
satisfy 

(1.31) d t p = —w 1 dip — pdjib 1 , 

(1.32) cW = -w I d I w J - -d J p{p) + g J , 

P 

(1.33) A$ = 4(p-A), 
where 



(1.34) § --^y T 3«^-4 9 $ -r 

A Newtonian potential can be defined by 

(1.35) $ = - + ^/jccV + —irprSux 1 ( p w J (fx. 

This potential satisfies the Poisson equation 

(1.36) A$ = p, 

while the acceleration due to gravity g J takes the familiar form 

(1.37) g J = -<9 J $. 

Together, equations (|1.31|) . (|1.32|) . (|1.36|) . and ()1.37j) show that solutions to the cosmological Poisson- 
Euler equations determine solutions to the standard Poisson-Euler equations on the covering space M. 

1.1. Notation. Before proceeding, we first fix our notation and introduce a number of function spaces 
that will be used in this article. Given a finite dimensional vector space V, we let H S {V) denote the 
standard Sobolev space of V- valued maps on T 3 . When V = M, we just write H s . The only two vector 
spaces that will be used in this article are M. N and the space of symmetric matrices Sn = { (u 1 -') G 
M NxN | «« = }. 

We denote the projection operator onto the L orthogonal complement of the constant function 1 by 

O(V0 =ip-(l\ip) L 2l V^eL 2 (T 3 ). 

Given {e a }^ =1 any basis for V, we use this projection to define 

N 



h s (v) = { v = ^ e HS{ y) | = o} 



a=l 



4n the Newton-Cartan theory, the fluid velocity 3- vector w 1 is the spatial part of a 4- vector w = dt + w T di [12]. The 
formula 11. 291 1 follows from the calculating the spatial components of w = ip*w. The other two formulas 11.281 1 and 111 .301 
follow from the definition of the pullback, i.e. p = ip*P an <i * = 
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We also define the standard hyperbolic evolution spaces 

£+1 

X TAs (V) = p| C*([0,T),H-*(V)). 

p=0 

and write Xx,e,s if V = K. 

1.2. Main result. The main results of this article is the following Theorem which establishes the 
existence of a wide class of 1-parameter families of solutions to the Einstein-Euler equations that converge 
in the limit e \ to solutions of the cosmological Poisson-Euler equations (|1.20|) - (|1.22|) . See Section [5] 
for the proof. 

Theorem 1.1. Suppose K > 0, n,sj € Z> , s > 3 + £, a £ H s , w' Q E H S (R 3 ), u 1 ^ € H S (S> 3 ) 
, u IJ S H S+1 (E>3), To as defined in Provosition 14.11 and o~q > is defined by (| 1 . 23[) where po — 
(4Kn(n + l))~ n aQ n . Then for eo small enough, there exists a T 6 (0,T ) independent of e € (0, eo), and 
maps 

Ue 3 € X T ^+i, s+ i(S 4 ) < e < e , 
a £ e X T>As , «;<(*) € X T/ , S (R 4 ) < e < e , 

,£+2,s+2, 

<reC u ((-e ,e ),R >tTo/2 ), a(0) = a , 

1 



y y e^((-e ,e ),M a ), y J (0)=-— / pou^x 



0o Jt 3 



such that 
(i) 



Q!e|t=0 = a 0, 



K)| t=0 = K,^), 



u°° \y J {e)+u° t J 
\y J (e)+uf eii IJ 

d t u° e J \ t =o 

t=o 



t=0 



where 



1„ „rn 3 



ft fi ™lt=o= - ^ ln(a)'^° - -a'<5 7J eii /J , 

^f| t=0 = - \ Ha)' (~y J (e) + , 

(ii) {u l ^> (t,x),a t (t,x),w\{t,x)} determines, via the formulas (jl.lip . (|1.13p . and (| 1 . 1 5|) . a 1-parameter 
family of unique solutions to the Einstein-Euler equations (jl.ip in the harmonic gauge on the 
common spacetime region (t,x) € M = [0, T) x T 3 , 

(hi) x), /5(i, x), w 1 (t, x)} solves the cosmological Euler-P oisson equations (|1.20|) - (|1.22|) on the 

spacetime region M , and 

(iv) there exists a constant C > independent of e € (0, eo) smc/i i/iai 

e^ll^W-lll^-x+ll^^-^WIU.-a 

+ ||Pe(t)-p(*)llu.-i <^e, 
- ^'(tJH^ + \\d K vii{t)\\ H s-i < Ce, 
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and 

\\ed t u™(t)\\ m -i <Ce, 
\\ed t u 0J -u° Q J + ^-u 0J \\ H s-i <Ce, 



2a 

''ii., // 

for all (t,e) E [0,T) x (0, e ), where 



ed t u IJ -uE\\ H s-i <Ce, 



Sg = 25Z8- 



2sM } -4c J 



2 

—t 
fi 

and fi, a, and ( J are defined by equations (|1.24|) . (|1.25|) . and (|1.26[) . 

2. The Einstein-Euler equations 

In this section, we adapt the formulation used in [10] to write the Einstein-Euler equations in a 
form suitable to analyze the limit e \ for initial data that does not satisfy (|1.2[) - (|1.3[) . The main 
change from the formalism used in [10] is that scale factor a e (t) appearing in the FLRW metric is now 
e-dependent through the function c(e). This additional freedom allows us to construct initial data that 
is more general than that in [10]. 

2.1. Reduced Einstein Equations. To derive a suitable symmetric hyperbolic system for the gravi- 
tational field equations, we introduce new coordinates related to old ones by the rescaling 

x° = x°/e, x J = x J , 

and let 

Q. = JL_ 

1 dx 1 ' 

In the new coordinates, the spacetime metric <?y and the FLRW metric hij (see (|1.4|l ) are given by 

Sa = JiJ-Qij and hij = + aS( 5 J % 5 U 

respectively, where Jf is defined in (|1.12p . The non-zero independent components of the Christoffel 
symbols 7k and the curvature H-ijki of the metric h^ are: 



(2.1) 7o/ 



e a'(t) 
2 a '■ 



(2.2) 7?z = ^a'(t), 

. - e 2 2aa" - (a') 2 

(2.3) Tioioi — —-, , 

4 a 

and 

e 2 

(2.4) 7^2121 = ^1313 = 7^2323 = ~^( a ') ■ 

As discussed in the introduction, we take the symmetric 2-tensor u v as our primary gravitational 
variable where 
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and 

(2.6) \g\ = -det(<7y) and \h\ = -det(hij) = a 3 - 

The metric can be recovered from the « iJ by the formula 

^-\h\det(g kl ) 

where 

(2.8) g ij = h ij + e 2 u ij . 

Substituting (|2.7p in to the standard formula for the Christoffel symbols gives 

(2.9) Y% = % + e 2 (-g l{l D 3) u kl + \g kl g im g jn DiU mn - \g k < %g mn Dfi mn + %9imS% i D j) Z lm ) , 

where (gij) = (g 13 ) and D k is the hij covariant derivative. Using this formula, the Einstein tensor 
G u of the metric gij is given by 

(2.10) \g\&' = € -\h\ [g M D k Diu ij + e 2 (4 + 4 + a$) + W + <<■( + e 2 ^ + 4e 2 T^] , 
where 

(2.11) a\ j = \(\g k i9mn - g km gin) (g rp g 3q - \g i3 g™)D p u kl D q u mn , 

(2.12) 4 = 2g kl (g n ^D m u^ l D n u km - \4 D rn u kn D n u ml — g mn D m u lk D n u jl ^j , 

(2.13) 4 = D k u ij DiU kl - D k v^ l Di%t? k , 

(2.14) b lj = g i] D k Diu kl - 2DiD k u k{i g j)l , 

(2.15) c l l = -(h ij eu kl + eu ir h kl )^n kl + ^U lk J l eu^ k h lm + ^K lk J*h^ k eu lm 

(2.16) 4 = -eu^eu kl \n kl + \ll lk J l eu^ k eu lm , 
and 

(2.17) T ij = (4^, C = ~/o- 
For later use, we define 

(2.18) 4 = -^iiujp—iisisl + ^n^sp + ^n km i {l h l)m 

so that 

ij ii -kl 
c l = c kl eu ■ 

To fix the gauge, we set 

(2.19) D t u tj = 0. 
For e > 0, it is clear from ()2.5|) that this is equivalent to 

D i (VW\S ij )=0, 

and this is easily seen to be equivalent to the harmonic coordinate condition 

fl ii (r?i-7?i) = o. 
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Defining the reduced Einstein tensor G 1 ^ by 

(2.20) G% =i - | (g kl D k D^ + e 2 (a« + a% + c#) + ecj j + e 2 4 3 + 4e 2 T«) , 
the Einstein equation G y = 2e 4 T ,y in the gauge (|2.19[) becomes 

(2.21) G%=2e 2 ^f l \ 

| ft. | 

where 

(2.22) f 4i = (p + e 2 p)v l v J + pg ij and u% = — ^. 

To write the reduced Einstein equations (|2.21|) in first order form, we introduce the variables 

(2.23) u ij = eu ij and v% = D k u ij . 
With these variables, we have that 

D k uf = D lUk 3 - -K kl J l u^ m , 

or equivalently 

(2-24) 8 k u? = ftujW + arSLuf™ - 27L"f m - f Wu**". 

In particular, this implies that 

(2.25) Bou'i = djv$ + 27<>f - 27& n «f n - ^o/ m ( V>\ 
and hence 

(2.26) = .9 00 a ^' + 2eu 0/ d / ^' + g 7 J a 7 w l j + edl 3 + e 2 4 , 
where 



(2-27) d^ = h kl ^X + ^L^ m ) 



and 



/n r>o\ j*7 fei I 1 -m in , ^ -U j)m\ . 01 f ^ - (i i)m ^ ~(i j)m 2 

(2.28) di=u I --7w< + -7L< )+« (-7j m «o ~ ~7om«r ~^ 



Using (|23 |l -([2T2 j) . we find that 



(2.29) = dg"*" W 



where 

a' / 3 



(2.30) 4 m =1 (-^« + $S$SjS? + h^flSTti - 
Setting 

(2.31) 47 = ^fC--7^ ) ^ 



and 

(2.32) /« = --^W^f, 
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equations ()2.25|) and ()2.26|) can be used to write the reduced Einstein equations (12.2111 in the following 
first order form 

-a(t)g 00 d t u l Q J = 2a(t)u 0I d I u% + -^a{t)g IJ d lU l ] + a{t)(c%u kl + 

(2.33) + ea(t) {a\ ] + c% + a% + 4 + <#') + 4ea(f) (V^ - |l|f « 

(2.34) a(t)g IJ d t v?j = ^(i)^^ + a(t)$"(e«7«* + 
and 

(2.35) = «« - *7&« i)fc - 

To be correctly defined, the reduced Einstein equations (|2.33[) - (|2.35[) require that the matrix c/ 4 -? is 
invertible. By assumption, <r(e) is analytic in a neighborhood of e = and cr(0) = <7o > (see (|1.23|l ). 
Thus there exists an eo > such tha10 

< y < (7(e) < 2a Ve G (-eo, eo). 



For fixed -^4/(3er ) < t < and t\ > 0, it is clear from (fl~9)) that 



:!/ '°ro + l) <a e {t) < (VtoEn + iy 



4 

for all (t,e) G [ro,ri] x (— eo,eo)- This implies that the set 

V = {(r«)|det(/i« + r«') > 0V(t,e) G [t ,ti] x (-e ,e )} 

is open and contains the origin (r u ) = 0, and moreover that the reduced Einstein equations (|2.33|) - (|2.35[) 
are well defined for all f G (tq,ti), e G (0, eo), and (eu v ) G V. We also note that 

2 f\g 



4e 2 ( f|T y - T iJ ) = 4(p - j*)d*c5g 



M 

(2.36) +e[8^W ) +^pft«« W ] + e 2 S J (e, ( 7, U A;i , /0 , W fc ), 

where the map Sq is analytic in all variables provided that (e, a) G (— e , eo) x (ao/2, 2<7o), and (eu y ) G V. 
2.2. Regularized Euler equations. In the coordinates (x l ), the Euler equations are given by 

(2.37) Vi? ij = 

where T lJ = (p + e 2 p)v l v^ + ptf J and the fluid velocity is normalized according to 

(2.38) s .jj* = — 1 

e z 

To derive a symmetric hyperbolic system for the Euler system, we follow the method of [1] and differ- 
entiate ()2.38|) to get 

(2.39) ViVjV* = 
which in turn implies 

(2.40) ffViVjtf = 0. 



The bound < < <r(e) < 2tro Ve 6 (— erj>eo) i s somewhat arbitrary and could be replaced by any bound of the 
form < ^S- < cr(e) < Cctq. However, as we are interested in the limit e \ 0, nothing is missed by assuming that C = 2. 
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Writing out (|2.37|) explicitly, we have 

(2.41) {d iP + e 2 B i p)v i v j + {p + e 2 p)(v j W l v i + + g lj B iP = . 
Next, we observe that the operator 

L\ = 5{ + e 2 v j Vi 

projects into subspace orthogonal to the fluid velocity v l , i.e. L\L\ = L 3 k and L\v % — 0. Applying this 
operator to project (|2.41[) into components parallel and orthogonal to v l yields, after using the relations 
{2~35|) - ([2T10]) . the following system 

(2.42) v%p +{p + e 2 p)L)ViV j = 0, 

(2.43) MijV k V k v j + —^L)d lP = 0, 

p + e A p J 



where 



Mij = gij + 2e 2 ViVj. 



(2.45) M^VkV 3 + ^-jLfaa = 



Changing to the Makino density variable a (see (jl.l5[l ) and multiplying (|2.42|) by the square of the 
function 

/= ( 1 + 4^TT)^ 2 

transforms the Euler equations (|2.42[) - (|2.43[) into 

(2.44) fv%a + ^//V,r' = 

a 

As discussed in the introduction, we need to introduce a new fluid four-vector by 

(2.46) w l = v l - C = v l - 

e 

So, letting 

(2.47) W M = (a,w l ) T 
allows us to write the system (|2.44l) - (|2.45l) as 

(2.48) A%d t W M = A^diWu + F M , 
where 

(2 49) A o_f/ 2 (l + -°) 

( } M ~{ 4%L? M«(i + e«,°) 



( 2 - 5 °) Al " = i-ku -m!J 



and 



(2 51) F M = 
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Next, a straightforward calculation using (12. 7p and (12.91) shows that 

Mij = hij + 2£°<5° + emij(e, a, u kl ,w k ), 
L)^S]-S^ + e£i(e,a,u kl ,w k ), 



3 o! 



It* 1 {^uf + 5 KL uf L ^ +e<&{e,oyiyi,w% 



Tki 



J kl v k v l = —Sjw 1 + 



-j J kl w k w l 



where the maps rriij, t\, qo, q±, and q J Q are analytic in all their variables provided that (e, a) G (— eo> eo) x 
(ctq/2, 2it) and (eu y ') G V. Using these expressions, we can decompose A M , A M , and Fm as 



(2.52) 
(2.53) 
(2.54) 

where 

(2.55) A M0 

(2.56) A Mfi 

(2.57) F M ,o 



Ki = Kifi + e4i(e, (T, U ij ,a, w% 

a m = a Ilo + <^M,i( e = CT ' uij > a > 

Fm = F M ,o + eFM,i(e, a, u lj ,u k \a, w l ), 



hij + 26?<5? 



2n u i 



» 3 

-(h ij +2S^)w I 



3a' 



(^ + 2<5?$) \-4° + i^°(3ug° - a^ L < L ) - (±u°° + 5 KL uf L ) - ^w 1 



and the maps A° M x , Aj^-g, -Fjw,i are analytic in all their variables provided that (e,cr) G (— £o7 £ o) x 
(cr /2,2cr ) and (etW) G V. 

2.3. A nonlocal symmetric hyperbolic formulation. To bring the reduced Einstein equations 
(|2.33[) - (|2.35[) into a form that is suitable to analyze the limit e \ 0, we replace the uj with the variables 

(2.58) Wf = u/ - 5 l 5 J di$ i: > 
where $ satisfies the Poisson equation 

(2.59) A$ = 4aII(p- fi). 

where A = 8 didj is the flat Laplacian. In addition to we also need the time derivative 

(2.60) $ = d t <f> 
which satisfies 

(2.61) A$ = 4all(d t (p - fj,)) + 4a'Tl(p - (j). 
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Using (|2.42l) - (12.43| to replace the time derivatives of p and dtw 1 in favor of spatial derivatives, we find 
that 



9t(p- M) 



3 a' 
2~a 



- M) - di(pw J ) + e 



ij 



(2.62) 



2 2 

e a 



Qo(e, c, u lj ,a, w l ) + Qi(e, a, u v ', a, w\ 9/a, diw 1 ) 



where Q v {y = 0, 1) are analytic in all variables for (e, er) e (— eo,£o) x (co/2, 2cto), and (eu^) £ V, and 
Qi is linear in (dia,diw l ). 
Letting, 



(2.63) 

we see from (12.62)) that (f> satisfies 

(2.64) 

Define 
(2.65) 



<t> = — <l|p — m)z 2 , 
e 



2 a 



L 2 ■ 



^M-UM, (m, dug), d^gD-(imi). nsn-iESi, szm-tirsn. ®m . nzm and dsn that w 

satisfies 



(2.66) 
where 
(2.67) 

(2.68) 

(2.69) 

(2.70) 

(2.71) 
(2.72) 

(2.73) 



A°d t W = -C 7 W + (Aq + eA[)diW + F + eFi, 



^40 _ / a g 



A 

G 





Am,q + e A°Ms) ' 

/a(l-eu 00 ) 0^ 

(<5 7J + eau IJ ) 

10 

\ 1/ 



A 1 



A 1 - 



0\ 



.4 





/2cm 07 au 7J 
au IJ 

\ 0/ 





1_ V A 7 ^ 

c 7 
c 7 



C G 

oy ' 

( 5 7J 0^ 

5 IJ 



\ 000/ 
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and 



(2.74) F = 



a(5i8i(u IJ didj$ + (/)- A P h H u kl ) - 85 { l pwrt + c 



u kl a m 



3 a / 
2 a 



<p|3^8°- 

3a 'a 



a 



5iju i q j ) 1 



2598?) 



(2.75) 4 



(2.76) Fi = 



For e > 0, 
gauge. It is 



( 



la' 
27 (P 



^°(3«8 



4na 



,00 



8klu? l ) - a ^8W) j 



a\a l { + a 2 " + a 3 J + 
±<S /j g>j($ - $ )<^ + au IJ [dj<S>5 l 5 J 
4(a 2 i|Q + Qi) 
V j(%,o — Fm,o) + Fm,i ) 

equation (|2.66|) is completely equivalent to the Einstcin-Euler equations in the harmonic 
this form of the Einstein-Euler equations that will be useful for analyzing the limit e \ 0. 



13 



2 



Jim M 
c Jkl u rn 



2.4. Well-posedness of the nonlocal system. The wcll-posedness of the non-local symmetric hy- 
perbolic system (I2.65P follows from the same arguments used in Section 2.4 of [10], and will not be 
repeated here. The well-posedness of the system (|2.65[) combined with its particular structure allows 
us to apply the local existence results of Schochet [13,14] (see also [3,4]) to obtain the existence of 
one-parameter families of solutions to (12.65|) on spacetime regions of the form M = [0, T) x T 3 where 
T is independent of e, and also to identify the cosmological Poisson-Euler equations (|1.20p - (|1.22p as the 
correct limit equations satisfied by the e \ limit of the solutions of (|2.65|) . The details of this are 
presented in Sections [4] and [5l 

3. Newtonian initial data 

In order to solve the initial value problem for the Einstein equations, we must first construct initial 
data that satisfies the following constraint equations on the initial hypersurface defined by t = 0: 



(3.1) 
(3.2) 

and 
(3.3) 



-2e 4 T 0j ")| 



= 0, 



A^| t=0 = o, 



gijv l v° 



0, 



4=0 



which are the gravitational constraints, the harmonic gauge condition, and the fluid 4-velocity normal- 
ization, respectively. 

To find 1-parameter families of solutions to these equations, we adapt the method developed by 
Lottermoser in [6]. We begin by writing (|3.2p as (see (|2.23|) ) 

kj I 



(3.4) 
which shows that 
(3.5) 

and hence, that 
(3.6) 



0. 



d.ju k 



kj I _ 



DiU 



kj I 

k lt=0 
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Substituting this into (12.14| yields 

(3.7) b°i\ t= =-g 00 ed t ui k . 
Setting 

H** = (g kl D k D lU v +6 ij ')| t=0 , 
it follows directly from (f2T26]> and ([3~7]) that 

(3.8) H oj = g IJ d lU °j j - g^edtu 3 / + 2evP I d I u° ] + ed[ j 
Setting j = J in the above expression while using (|3.4p gives 

(3.9) H 0J = g KL d K u^ ~ g^edtuj 1 - 2eu 0I d I uf + ed 
From the non- vanishing Christoffel symbols (|2.1|l - l|2.2|l . we get that 

(3.10) u° °=ed t u 00 , 

(3.U) - 0J - ^jrflJ i a ' „,0J 

(3.12) 

and 
(3.13) 



Uq =ed t u» J + — u UJ , 
2a 

' _ -J J , a ' 



Jtu H u 

a 



Oku -t a u , 



(3.14) 
(3.15) 



0J = d K u 0J 



' K 



i IJ 



;J „,oo 



— n, KJ -L- —y 7/ U 

2 U + 2a KU ' 



Using (O, and ([3~TU|) - ([3~15|) . a straightforward calculation (see P7l3|) and (|2~2"7|) ) shows 

that 



(3.16) 



e<9 4 w} J + e« + c UJ ) = d^tdtu 1 - 1 ) + e 



3 /a' 



2 V a 



,J0 



2a 



2a 2 



rdju 



01) 



Setting j = in (Ml) while using ([2~23]) and gives 



(3.17) 



2_(0 7)m 



t5 (0 -T)m 



2 U 0/ a 7U f + < u + 



(HI 



,110 



Next, we decompose the gravitational variables as follows 



(3.18) 
(3.19) 
(3.20) 



U°3 = lylSj+U ^, 



U IJ = eu IJ , 



d t u 



■ "0 > 



where we assume that the y are constants and 



(3.21) 



f j d 3 x = 0. 
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In terms of these variables, the harmonic conditions become 

1 



(3.22) 
(3.23) 



ed t u m = --d!U 10 
e 



d t u 0J = -d I u IJ -^ln(a)'(^y J 



-OJ 



Using (pJ0]l . (I3~9l) . ([3~T6]) . (f3TTT|> . (|3~22|) . and ([3~23]) . the gravitational constraints (|3~Tjl take the form 
(3.24) Au 00 - Aa{p- p) + e[f°(e,a,y,u, p,w) + f°(e,a,y,u,d%u,d x u,u ) + /° (e, a, y, u, d x u, u )] = 0, 

3 /aV 

,2 



OJ 



ea 



-d J u 00 



:(£L)V-W 



(3.25) 



/i J (e, c, y, P, + //(e, a, y, u, d 2 x u, d x u, u ) + / 2 J (e, a, y, u, d x u, u Q ) 



= 0, 



where (i) for any R > there exists an eo > such that the maps f a (a = 2, 3, 4) are analytic in all their 
variables provided |e| < e , er € (<J /2, 2cr ) and \u tJ \ < i?, (ii) /| is linear in (OkOlu 13 ^ku 13 ,d t u 13 ), 
and (iii) j 3 2 is quadratic in (<9 t u y , Dku 13 ). Here, a, a', and /i are determined by evaluating the formulas 
HU), III]) and (fTTUI) at < = 0, that is 

= 1, u = a, and a' = \ —a. 

V 3 



(3.26) 

Also, a calculation using (|2.7[) shows that (|3.3[) can be written as 
(3.27) w -ef (e,(T,y,u ij ,xv I ) = 0, 

where the map /o is analytic provided |e| < eo, o S (o"o/2, 2<7o), | z2' 1 - 7 | < i?, and |ur| < i?. 

, ^ J € ff s (S 3 ) , u IJ € Bfl(iJ s+1 (§ 3 )). 



Theorem 3.1. Suppose s > 3/2, R > 0, a € H s , w 1 £ H 

Then there exists an eo > 0, and analytic maps 

(-e ,e ) 3 e i 
(-e ,e ) 3 e 
(-e ,e ) 3 e i 

and 



a(e) e M> , 
/(e) S M 3 , 



-e ,e ) 3 e i 



= it 



f e /i s (l 



sucft i/iaf /or euen/ e 6 (— eo, eo) 



a|t=o = a, 
(wi)\ t=0 = (w° e ,w I ), 



where 



-,00 



V(e)+«. 



.70 



?/(e) + zI ( 



OJ 



-,/./ 



001 

« lt=o 



-9 7 ii 70 -^ln(a)'^ u --a' ( S 7J ezi iJ , 



OJ 



7/ /J 

"0 



1 o 
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and 



d t u° t J \ 



-diu 



i J 



>(a)' V(e 



-.OJ 



satisfy the gravitational constraint equations (|3.ip . the harmonic gauge condition (|3.2|1 . and the fluid 
velocity normalization (|3.3p . Moreover, the maps cr(e), w®, y J (e) and vPJ admit the expansions 



a(e) =a + O(e), w° e = 0(e), y J (e) = yi + 0(e), 
+ O(e), and w° J = 0(e), 



where 



and 



(TO 



T 3 



pd 3 x, 2/0 



CO Jt 3 



= A Hp-o-o)- 

Proof. Applying the projection operators II and I — II to the equations (|3.24p - (13.25p . while observing 
the averaging conditions (|3 . 2 1 [) and the definitions (|3 . 26[) . yields the equations 



o-- (l|p>L 2 + -|(l|-F (e, o-,y, u, u ,a, w)) L 2 = 0, 

y J + -{l\pw J ) L 2 - ^-e{l\F J (e,a,y 7 u,u ,a,w)) L 2 = 0, 
Am 00 - 4H(p) + eli(F°(e, a, y, u, u , a, w)) = 0, 



(3.28) 

(3.29) 
(3.30) 

and 

(3.31) Aw 0J + ell(-4pw j + eF J (e, a, y, u, u , a, w)) = 0. 

These equations together with (|3.27p form the complete set of constraint equations to be solved 

Using the same arguments as in Section 3 of [10], it can be shown that for any s > 3/2, R > 0, and 
Co > that there exists an e > such that the maps 

{e,a,y J ,u ij ,ul J ,a,w j ) 3 (-e ,e ) x {a /2,2cr ) x K 3 x 

B R (ti s+1 (§ 4 )) x H S (S 3 ) x H s x H S {R 4 ) h-> F' J € H s ~ l {^) 

and 

(e,a,y,^'V) 9 (-e , e ) x (a /2,2a ) x M 3 x B R (H S+1 (S 4 )) x tf s (M 3 ) ^ f e H s 
are analytic. Setting 

•0 = (u /J , Mg' 7 , a, w J ) T and 77 = (cr, y J ', u - 1 , u; ), 
allows us to write the constraint equations (|3.27p . (I3.28[) - (|3.31|) as the following single equation: 

(3.32) L(<i/j,n) + eM(e,il},r]) = 
where 



(3.33) 



( a-(l\p) L 2 \ 
y J + ±(l\pw J ) L 2 
Au°i - m(p)8 0: > 



\ 



J 
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and M(e,ip,rf) are analytic maps. It is clear that 

/ 

(3.34) 77WO - 



4<5^'A- 1 II( ( o) 
V o 



solves 
(3.35) 



LftMty)) = 0. 



Also, a straightforward calculation shows that 
(3.36) D 2 L(i(),ri(ip)) ■ 8r) 



I 5 ° , \ 

V Sw° J 



and it follows easily from this equation and the invertibility of the Laplacian A : H s+1 
the linear map 



H" 



that 



(3.37) 



D 2 L(ijj, rj(ip)) :Rx 



x H 



s+l. 



x H s 



x H S -\R 4 ) x H s 



is an isomorphism. The proof of the Theorem now follows from (|3.32[) - (|3.37p and the analytic version 
of the Implicit Function Theorem (see [2], Theorem 15.3). □ 

4. Limit equations 

The evolution equations (|2.65|) are now in a form to which the theory of singular symmetric hyperbolic 
equations [4, 13, 14] pioneered and developed by Kreiss, Klainerman, Majda, and Schochet applies. We 
know from the theory developed in these works that the appropriate limit equation that is satisfied by 
solutions of (|2.65[) in the limit e \ is 

(4.1) A°d t W = A^djW + F + c 7 w 

where 



Fr> — Fn 



A 



e=0' 







^o : a/| c =o 



and 



(I 



A 



'a 0\ 
S IJ 
10 
\0 1/ 

As we shall shortly see, solutions to the cosmological Poisson-Euler-Makino equations p.l6[) - (j!.18p 
determine solutions to the limit equations (|4.ip . But first, we record the following local existence and 
uniqueness result for the Poisson-Euler-Makino equations. 

Proposition 4.1. Suppose s > 3 + £, &q E H k , and Wq £ H k (M. 3 ). Then there exists a maximal time 
To and a unique solution 

a e X To ,i, s , w 1 € X ToAs (R 3 ), $ e X To . e+2 . s+ 2 
to the Poisson-Euler-Makino equations (|1.16|) - (jl.l8|) with initial data a\t=o = <5q an d w 7 |t=o = ^o- 
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Proof. Since the Poisson-Euler-Makino equations (|1.16|) - ljl.l8|) form a (non-local) symmetric hyperbolic 
system, the proof follows from standard theory. For example, see [15], Chapter 16. □ 

Proposition 4.2. Suppose {a, w 1 , $} is the solution to the Poisson-Euler-Makino equations from Propo- 
sition |4~T) and let 

a = exp ( f (— f p(s)d 3 x J ds \ , 



The 



T 3 



ft 


= / pd 3 x, 






= / pw J d 3 x. 




~ij 


= 2*4' (- ( 


¥) 






4 


wf 


= 5u 0/ , 






= 0, 






a 

= 0, 











2d p^ J 



ha = -S?$ + a8!6?6 I j, 

4f - 1 (-f 4« + WW + ^W^*? - WW) > 



W=(4 j ,W i f i ,u ij ,4>,a,w i ) T , 
w = (u$,J/, 0,0,0, Of, 

defines a solution to the limit equations (14. lj) on t/ie spacetime region Mq = [0, To) x T . 

Proof. The proof follows by a straightforward calculation that verifies W is a solution to the limit 
equation (|4.ip . □ 



5. The Newtonian limit 
We are now ready to prove Theorem ll.il 
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Proof of Theorem\Tj} Given <5 G *o € iJ s (R 3 ), G H s {§ 3 ) , u IJ G iP +1 (§ 3 ), we let 

a e \t = o = cto, 
K)|t=o = K>o). 

qj \ / € 



cm 



<9 t u 



t=o 



t=0 U 

3 , , ,/„no 



where 



r\ -< 

o t u 



1 



- dl u{°- l\n(a)'u ,J e tJ --a'6 I jeu 1 \ 



-,1-J 



e \t=0 



d ^e J \ t=0 = ~ d rt" - \ ln(o)' ( V(e) + fi 



OJ 



be the initial data from Theorem 13.11 By construction, this data solves the constraint equations (|3.1[1 - 
and depends analytically on e, and satisfies 



a c = (To + 0(e), u£ = 0(e), /(e) = y J + 0(e), 
> + 0(e), and u° J = 0(e), 



where 



and 



pd 3 x, yo 



= a Hp- c ) 

A straightforward calculation then shows that 
(5.1) 
satisfies 



(5.2) 
where 



W e {t) = (u l l(t),Wrjt),unt),Mt)^e(t),wi(t)Y 

W \ t=0 = (ul j U=o, 0, u ij \ t=Q , 0, oo, 



= (4fe(n + l))-V n , 
a = exp [ / ( — f p(s)d 3 x I e?s 



T 3 



pd 3 x, 



T' 3 



T 3 



«y-*M(-(±c)'-^) 
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and {a,!!! 1 , $} is the solution to the Cosmological Poisson-Euler equations from Proposition 14. II that is 
generated from the initial data <5|t=o = So and w / |t = o = w^. Also, it can be verified using the evolution 
equations (|2.65[) that our choice of initial data guarantees that 

(5.3) ||d t W 6 | t=0 || H . <1 0< e <e 

where 

H s = (tf s (§ 4 )) 3 x E x H s x H S (R 4 ). 

Together, Proposition 14.21 and (|5.2|) - (|5.3() allow us to apply Theorem 1 of [14] (see also the remarks in 
Section 1 of [14]) and conclude (shrinking eq if necessary) that for any T < To, there exists maps 

£+3 

W t S p| C p ({0,T),H s ~ p ) 0<e<e o , 

p=0 

such that 

(i) W e (t, x 1 ) solves equation (|2.66p on the spacetime region (t = x°, x 1 ) € D = [0, T) x T 3 , 

(ii) W e satisfy the estimates 

\\We(t)\\w + \\dtWe(t)\\ H .-i <1, 

for all (t,e) £ [0.T) x (0.e ), and 

(iv) 

\\W e (t)-W(t)\\ H .-i<e, 
for all (t,e) € [0,T) x (0, eo), where W is the solution of the limit equations from Proposition 

EH 

Finally, we observe that similar arguments used in the proof of Proposition 6.1 in [8] show that 

determines, via formulas (| 1 . 1 5[) and Ql.ll[l - (|1.13[l . a solution to the Einstein-Euler equations (|1.1[) in the 
harmonic gauge (|2.19[) . and moreover, that 

ed t uf 

ed t u 0J 

U 



00 




'0,ei 




0J 


a -0J 


0,e 


2a e 


/./ 


a ' -IJ 


0,e 




a 



ed t u 
and 

d K u™ = W K <e -a'eu° e K + d K ® e , 
d K u° J = W% e - a -eu? J + ^6' K eu™, 

<JKU e — W Ke — —ed K u e ' . 
This combined with the statements (i)-(iii) above completes the proof. □ 
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